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The total mass and dark halo properties of the Small 
Magellanic Cloud 
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ABSTRACT 

We discuss the total mass and dark halo properties of the Small Magellanic Cloud 
(SMC) for reasonable F-band stellar-to-mass-to-light ratios (Ms/Ly) based on the 
high-resolution neutral hydrogen (HI) observations of the SMC. We adopt both the 
Burkert and the NFW profiles for the dark matter halo of the SMC in order to model 
the radial profile of the observed rotation curve. We show that the mass (Mdm) and 
the mean density (pdm) of the dark matter halo within the optical radius (^ 3 kpc) 
of the SMC can be significantly low for AIs/Ly = 2 — 4 reasonable for a galaxy 
dominated by older stellar populations. For example, the estimated Mdm and pdm are 
7.5 X 10'^Mq and 6.7 x IQ-^'Mq pc'^, respectively, for M^/Ly = 3.8. The models 
with lower Mg/Ly (< 1.0) can show higher Mdm and pdm yet have difficulties in 
reproducing the inner rotation curve profile both for the Burkert and the NFW profiles. 
The Burkert profile with a larger core radius (> 1 kpc) (thus a low density) and a 
large mass (Mdm > 3 x IO^Mq) can be better fit to the observed slowly rising rotation 
curve than the NFW one for a reasonable M^/Ly ~ 2. This means that most of the 
dark halo mass can be initially located outside the optical radius in the SMC and 
thus that the dark halo would have already lost a significant fraction of its original 
mass due to the strong tidal interactions with the Galaxy and possibly with the Large 
Magellanic Cloud (LMC). We suggest that the dark mater halo of the SMC is likely 
to have the initial total mass and core radius as large as, or larger than, 6.5 x IO^Mq 
and 3.2 kpc, respectively. We discuss limitations of the present model in estimating 
the total mass of the SMC. 

Key words: galaxies: dwarf - galaxiesistructure - galaxiesikinematics and dynamics 
- galaxies: star clusters 



1 INTRODUCTION 

Observational studies of the neutral hydrogen (HI) gas in 
the Magellanic Clouds (MCs) have played a key role in the 
understanding of many aspects of their formation, such as 
the formation of the Magellanic Stream (MS; Mathewson 
et al. 1974; Putman et al. 1998), gaseous kinematics of the 
Large Magellanic Cloud (LMC; Staveley-Smith et al. 2003), 
and the formation of the Magellanic Bridge (MB; MuUer 
et al. 2004). These observations furthermore have attracted 
much attention from theoretical studies of the roles of LMC- 
Small Magellanic Cloud (SMC)-Galaxy interactions in the 
formation of the MS (e.g., Murai & Fujimoto 1980), the 
importance of the last strong LMC-SMC interaction in the 
formation of the pecuhar HI morphology (Bekki & Chiba 
2007), and the origin of the bifurcated kinematics in the 
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MB (e.g., Muller & Bekki 2007). Recent observations of the 
structural and kinematical properties of the HI gas in the 
SMC have provided valuable information on the mass dis- 
tribution of the SMC (Stanimirovic et al. 2004, S04). 

One of intriguing results by S04 is that the HI gas in the 
SMC has rotational kinematics with the maximum circular 
velocity (14) of ~ 60 km s~^. This kinematics is in a striking 
contrast to the stellar kinematics with no or little rotation 
(e.g., Harris & Zaritsky 2006), which recently have led Bekki 
& Chiba (2008) to suggest that the SMC have experienced a 
major merger event long time ago. Another intriguing result 
in S04 is that 14 ~ 60 km s~^, corresponding to the total 
mass of 2.4 x 1O''M0, can be explained by the mass distri- 
bution of baryonic components (i.g., gas and stars) without 
a dark matter halo for the central 3 kpc: this result however 
depends strongly on the adopted stellar distribution. These 
properties certainly make the SMC a highly unusual dwarf 
galaxy. However, no detailed dynamical studies have been 
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Table 1. Parameter vales for the representative models. 

Md„, (M0)£ ro(kpc)^ 



Vlodel_^ 
NDMr 


40 


A/s/Lv 
3.8 


NDM2 


30 


3.8 


NDM3 


40 


2.3 


Bl 


40 


2.3 


B2 


40 


2.3 


B3 


40 


0.6 


B4 


40 


1.0 


NFWl 


40 


2.3 


NFW2 


40 


2.3 


NFW3 


40 


2.3 



6.5 
79.3 
79.3 
79.3 
6.5 
6.5 
79.3 



3.2 
9.5 
9.5 
9.5 



(kpc) /. 



2.5 
5.1 
17.7 



10 

5 
5 



Comments 

without dark matter halo 

lower inclination angle 

lower A/s/Lv 

Burkcrt profile 

very low Afs/Lv 

NEW profile 



^ The mass model for the dark matter halo of the SMC: "NDM", "B", and "NEW" represent a model with no dark matter (or with no specific 
dark matter profiles considered), the Burkert profile for the radial density profile of the dark matter halo, and NEW one, respectively. 

The inclination angles of the gas disk in the SMC. 
'^ The V-band stellar-to-mass-to-light-ratio. 
'^ The total mass of the dark matter halo in units of 10^ Mq. 
^ The scale radius in the Burkert profile. 
J The scale radius in the NEW profile. 
9 c {= rvir/^s) parameter in the NEW profile. 



conducted to explain the origin of such unusual dark-matter 
halo. 

The purpose of this paper is thus to discuss the exis- 
tence of dark matter within the central 3 kpc of the SMC, 
based on the detailed models for the mass distributions of 
the dark matter halo and the baryonic components. We de- 
rive the possible total mass (Afdm) and mean density (pdm) 
of the dark matter halo within the optical radius of the SMC 
required to reproduce the observed HI rotation curve. We 
discuss whether the modeled distributions of the dark mater 
halo in the SMC can be consistent with the observed slowly 
rising rotation curve. We also discuss a possible initial mass 
of the SMC before the interaction with the Galaxy and the 
LMC. 

This paper focuses on implications and interpretation 
of observational results presented in S04. We stress that 
the observed rotation curve of the SMC depends on sev- 
eral adopted parameters: inclination, position angle, the disk 
thickness and shape of the SMC. In addition, possible radial 
gas- flows in the SMC, caused by the LMC-SMC-Galaxy in- 
teractions, can result in non-circular motions of HI gas and 
make the estimate of the rotation curve even harder. Some 
of the adopted parameters can not be directly determined 
by observations, because the 3D HI structure remains ob- 
servationally unclear for the SMC. Since observational un- 
certainties involved in the derivation of the rotation curve 
have already been discussed in S04, we just briefly summa- 
rize them in this paper and explain limitations they impose 
on the derivation of the dark matter halo. We also point 
out that the total mass of the SMC derived in S04 agrees 
well with the estimate based purely on the stellar kinematics 
(Harris & Zaritsky 2006). 



2 POSSIBLE TOTAL MASS AND MEAN 

DENSITY OF THE DARK MATTER HALO 
FROM OBSERVATIONS 

We first derive the possible mass and density of the dark 
matter halo of the SMC from the observed rotation curve 



(S04) without using models for the dark matter in §2.1 and 
§2.2. We then try to find the density profile of the halo 
that can fit well to the observed rotation curve using models 
predicted from previous theoretical studies in §3. We discuss 
how uncertainties in the observed HI rotation curve (e.g. a 
possible existence of radial gas infiow) affect our results in 
Section §2.3. 



2.1 The observed rotation curve 

Before discussing the possible total mass and density of the 
dark matter halo of the SMC, we firstly describe the ob- 
served radial profiles for the total mass (M) and the rota- 
tional velocity (K) of each component (e.g., gas and stars) 
within a distance (R) from the center of the SMC based on 
the results by S04. We here assume that the SMC has a 
stellar luminosity of 4.3 x lO^L© (i.e.. My — —16.76 mag), 
which is significantly smaller than that in S04: the rotation 
curve from stars in the present paper is different from that 
in S04. 

As described in S04, the tilted-ring algorithm was used 
to derive the HI rotation curve. Keeping the kinematic cen- 
ter fixed, and assuming an inclination angle of 40 degrees, 
the final rotation curves for the receding and approaching 
sides of the HI velocity field were derived. As the line-of- 
sight velocity dispersion in the SMC is high (mean value 
22 km s~^), the asymmetric drift correction was applied, 
resulting in the final rotation curve being higher by about 
10 km/sec for R > 0.5 kpc than the observed one. In addi- 
tion, to derive the deprojected rotation curve due to the HI 
gas alone, an exponential density law was assumed for the 
vertical HI disk distribution with a scale height of 1 kpc. 
The deprojected rotation curve due to the stellar compo- 
nent was derived from the surface brightness distribution at 
V-band. For full details please see S04. Thus the circular 
velocity (K) at a given radius is a combination from dark 
matter (Vcdm), stars (Vc,s), and gas (Vc,g). Fig. 1 shows de- 
rived radial profiles of M and K for the model with the 
most plausible value for inclination {i — 40), M^/Lv = 3.8, 
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and with no dark matter (referred to as "NDM" model from 
now on). One source of uncertainty that has not been taken 
into account previously is the inclination angle. S04 esti- 
mated i at each radius and found that the mean value of 
i is 40° ± 6° within the central 2.3 kpc. After this radius i 
was found to rise slightly outward. We have hence assumed 
the observational error (Ai) in i (at a given R) A = 5° and 
propagated this value to estimate the uncertainty in M and 
Vc (as shown in Fig. 1). Later, we investigate the models 
with different mean i to infer the minimum and maximum 
possible masses of the dark matter halo of the SMC within 
the central 3 kpc. 

It is clear that the observed profiles are quite similar to 
the profiles for the total component (i.e., the combination of 
stars and gas) , though M and K around R — 1 -- 2 kpc for 
the total component are slightly smaller than the observed 
ones. Clearly, the observations can be well reproduced with- 
out invoking the dark matter halo for the optical radius of 
the SMC. There is only a slight difference between the total 
mass (Mt) and the dynamical mass estimated from K within 
the central 3 kpc. If this difference is due to the presence of 
the dark matter halo, the total mass (Mdm) and the mean 
density (pdm) within the central 3 kpc are 7.5 x lO^M© and 
6.7 X 1O"*M0 pc"^, respectively. 

What would be a reasonable value for Ms/Ly based on 
integrated optical colors? We use the formular described by 
Bell & de Jong (2001, BDOl) and the observed colors of the 
SMC: B-V (=0.41) a.ndB-R (=0.70). If we use the low 
metallicity models with Z = 0.008 in BDOl, the most plau- 
sible Ms/ Ly is about 1. However, M^/Ly can range from 0.6 
to 1.0 for different stellar population models and star forma- 
tion histories (BDOl, S04). The stellar population synthesis 
models for [Fe/H] = -0.68 and ages of 2, 10, and 13 Gyr by 
Vazdekis et al. (1996) show M^/Lv = 1.1, 3.8, and 4.4 re- 
spectively. This means that Ms/Ly strongly depends on the 
luminosity-weighted age of the stellar populations. As it is 
observationally not clear what are the dominant ages of stel- 
lar populations, we assume that Ms/Ly is a free parameter 
in the present study. Given that the SMC contain a signifi- 
cant fraction of older stellar populations (Harris & Zaritsky 
2004), we consider that the reasonable Ms/Lv would be 2 — 4 
for the observed large fraction of old stars in the SMC (e.g., 
Harris & Zaritsky 2004). 

2.2 Uncertainties in Ms/Ly and i 

We now explore how Ms/Ly and i affect Mdm and pdm- 
Generally, for a given Ms/Ly, i < 40 degrees results in sig- 
nificantly larger Mdm and pdm- We first calculate Mdm from 
the total mass (derived from the observed Vc at _R = 3 kpc 
for a given i) and and the observed gaseous (Mg) and stel- 
lar masses (Ms) for a given Ms/Ly. Mdm = M — Mg — Ms 
so that no cold and warm gas possibly in the SMC's halo 
can be considered. We then estimate pdm from the derived 
Mdm within the central 3 kpc of the SMC. Fig. 2 shows 
that Mdm and pdm are significantly smaller than Mb (i.e., 
the total mass of the baryonic components including both 
gas and stars) and pb (the mean density of the baryonic 
components), respectively, for a reasonable i (= 40°) and 
1.5 ^ Ms/Ly ^ 3.8. Only if i is largely underestimated 
by ^ 10° in the observation by S04, Mdm and pdm can be 
significantly large: the model NDM2 with i — 30° shows 




R (kpc) 

Figure 1. The total masses (M, upper) within R (where R is 
the distance from the center of the SMC) and circular veloc- 
ities (Vc, lower) at R for all components (solid), dark matter 
(DM) one (dotted), stellar one (short-dashed), and gaseous one 
(long-dashed) in the model NDMl without dark matter. Here the 
inclination angle i and the V-band stellar-to-mass-to-light-ratio 
(Ms/Ly) are 40° and 3.8, respectively. For comparison, the ob- 
served values arc plotted by filled circles. 



Mdm = 1-6 X IO^'Mq and pdm = O.OISMq pc~^. Since i is 
determined by the tilted ring algorithm at each R in S04, it 
is unlikely that i is underestimated by ~ 10° for the entire 
region of the SMC. We thus conclude that the larger Mdm 
due to the smaller i is highly unlikely. 

It should be stressed here that (i) the derived total mass 
of 2.4x lO^M© for the best i (= 40°) in S04 is consistent with 
that derived from stellar kinematics of the SMC (e.g., Harris 
& Zaritsky 2006) and (ii) the best i is also consistent with 
that for the best theoretical model that can explain phys- 
ical properties of the Magellanic Stream and Bridge (e.g., 
Yoshizawa & Noguchi 2003). These strongly suggest that al- 
though i and position angles change with R in S04, the best 
i and position angle in S04 can be very close to the true 
ones. This implies that the adopted tilted-ring algorithm in 
S04 is still useful in accurately estimating the rotation curve 
for the inner part of the SMC (i? < 3 kpc), where the infiu- 
ences of the LMC-SMC-Galaxy interaction can be relatively 
weak. 

The model NDM3 with Ms/Ly = 2.3 shows Mdm = 
7.8 x IO^Mq and pdm = 6.9 x 1O"^M0 pc"^. Mdm and 
pdm are smaller for a larger Ms/Ly. Mdm and pdm are 
4.3 X IO^Mq and 3.8 x lO'^M© pc"^, respectively, for 
Ms/Ly = 3.1. The total mass of the dark matter halo within 
the optical radius of the SMC in the model NDM3 is, how- 
ever, still low, given that the total baryonic mass of the 
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Figure 2. The mean density (pdmi upper) and the total mass 
(Afdmi lower) of the dark matter halos within the central 3 kpc 
of the SMC as a function of Ms/Lv, for i = 40° (solid), i = 30° 
(dotted), and i = 50° (dashed). For comparison, the mean density 
(Pdm) a-ud the total mass (M^,) for the baryonic component of 
the SMC for Mg/Ly = 3.8 arc shown by dashed lines. This figure 
indicates the possible upper and lower masses of the SMC dark 
matter halo. 



SMC is 1.6 X lO^M©. The SMC may well be a baryon- 
dominated dwarf galaxy. These results imply that only if 
the SMC is dominated by a very young stellar population, 
and thus has a small Ms/Ly, the SMC can have a certain 
amount (~ lO^Mo) of dark matter within its optical radius. 
As a conclusion, the SMC can contain a certain amount 
of dark matter within its optical radius, if it has a signifi- 
cant fraction of young stellar populations and thus a small 
Ms/Ly (~ 0.6): Mdm and pdm depend strongly on Ms/Ly 
such that they are both larger for smaller Ms/Ly for a given 
i. The SMC can hardly contain dark matter, if it is dom- 
inated by old stellar populations with ages larger than 10 
Gyr and thus has Ms/Ly ^ 4. In the next section, we dis- 
cuss possible mass distributions for the SMC's dark matter 
halo that can be fitted to the observed K profile by assuming 
Ms/Ly = 2.3 (corresponding to ages of stellar populations 
similar to 5 Gyr) that would be reasonable for a galaxy hav- 
ing a large fraction of older stellar populations like the SMC. 



survey of 2046 red giant stars has suggested that the older 
stellar components of the SMC have a velocity dispersion 
(a) of --^ 27.5 km s~^ and a maximum possible rotation of 
~ 17 km s^^ (Harris & Zaritsky 2006). The older stars are 
observed to have a more regular distribution and appear 
to be a slightly flattened ellipsoid rather than a disk (e.g., 
Cioni et al. 2000). These implies that the SMC has a stellar 
spheroid rather than a rotating disk: the SMC is unlikely to 
have a bar composed of older stars. The lack of a stellar bar 
would not cause the non-circular motion of gas in the SMC 
and thus suggests that the observed rotation curve by S04 is 
highly unlikely to be influenced by the non-circular motion 
of gas due to the bar. 

As discussed by S04, numerous expanding shells in the 
SMC can provide some non-circular motions and thus be the 
main reason for the possible non-perfect rotation curve of 
the SMC. S04 however considered this effect (i.e., introduc- 
tion of velocity dispersion) by calculating the asymmetric 
drift correction. The last LMC-SMC tidal interaction about 
0.2 Gyr can give tidal perturbation to the SMC and thus in- 
fluence gas dynamics of the SMC in its outer part (Bekki & 
Chiba 2007). Although the outer part of gas (3 < i? < 5 kpc) 
can be stripped to form the Magellanic Bridge during the in- 
teraction (e.g., Yoshizawa & Noguchi 2003; MuUer & Bekki 
2007), the inner part of the present SMC (i.e., 0.2 Gyr after 
the last interaction) can still show clear rotation. Therefore, 
if the disordered gas motion due to the last LMC-SMC in- 
teraction can influence the present mass estimation of the 
SMC from the observed rotation curve, the influence would 
be signiflcantly small: a possible influence would be slight 
underestimation of the mass due to Vc smaller than the true 
value resulted from a higher velocity dispersion (caused by 
the last interaction) in the outskirt of the SMC. 

The 3D structure (e.g., thickness and shape) of the HI 
gas in the SMC remains observationally unclear. We thus 
can not make robust conclusions on how the thickness and 
shape of the HI can influence the mass estimation of the 
SMC. S04 adopted a reasonable assumption on a constant 
vertical scale-height (Ikpc) and considered asymmetric drift 
corrections in deriving the rotation curve of the SMC. We 
here suggest that a possible large intrinsic ellipticity (e.g., 
e ~ 0.3) of the HI distribution can introduce an observa- 
tional error in the estimation of the rotation curve in the 
SMC. Since the possible range of e for the SMC is observa- 
tionally unclear, it is hard for the present study to make a 
quantitative estimation for the possible observational error. 
Keeping these limitations in mind, we discuss the models for 
the dark matter halo of the SMC in the following sections. 



THE MODELS FOR THE DARK MATTER 
HALO 



2.3 Limitations of the mass estimation 

We have so far considered that the observed rotation curve 
is very close to the true one with some possible uncertainty 
of the disk inclination (j). It should be here stressed that 
the observed rotation curve can be slightly different from 
the true one owing to the possible non-circular flow of gas 
caused by the presence of a stellar bar in the SMC. The latest 



3.1 The NFW profile 

We now investigate the M- and K-profiles in the models 
with dark matter halos and Ms/Ly = 2.3 reasonable for the 
SMC having a large fraction of intermediate-age and old stel- 
lar populations. We adopt the same radial density profiles 
(described later) of dark matter halos as widely adopted in 
other studies (on isolated systems with no apparent galaxy 
interaction) , though the SMC appears to be interacting with 
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Figure 3. The relation between the core radius (ro) of the SMC's 
dark matter halo described as the Burkert profile (B95) and the 
total mass of the SMC (Mt) for Ms/Ly = 2.3 (solid), Ms/Ly = 
3.1 (dotted), and Ms/Lv = 3.8 (dashed). 





R (kpc) 

Figure 4. The same as Fig. 1 but for the model Bl in which 
the Burkert profile with ro = 3.2 kpc, i = 40°, and Mg/Ly = 2.3 
are adopted. Mdm and Mt are 6.5 X lO^M© and 8.1 X IO^Mq, 
respectively, in this model. Note that the Burkert profile with a 
large core radius can well reproduce the slowly rising Vc profile. 



the LMC and the Galaxy. We think that this adoption can be 
reasonable, because we discuss the density profile in the in- 
ner part of the dark halo, where the influence of the Galaxy- 
LMC-SMC interaction is relatively small. 

We adopt the following two profiles for the radial den- 
sity distribution of the dark matter halo of the SMC: the 
"NFW" profile with a central cusp predicted from hierar- 
chical clustering scenarios based on the cold dark matter 




R (kpc) 

Figure 5. The same as Fig. 1 but for the model NFWl in which 
the NFW profile with r^ = 2.5 kpc and c = 10, i = 40°, and 
Ms/Ly = 2.3 arc adopted. M^^ and Mt are exactly the same as 
those used in the Burkert one shown in Fig. 4. Note that the fit 
of the NFW profile to the observation is worse than the Burkert 
profile in Fig. 4, though the NFW profile can better reproduce 
the Vc profile in the outer part (R ~ 3 kpc). 



models (Navarro et al. 1995) and the "Burkert" profile with 
a flat "core" (Burkert 1995; B95). 
The NFW profile is described as: 



Pdm('") = 



PO 



(r/r,)(l + r/r,)2 



(1) 



where r, po, and rs are the distance from the center of the 
cluster, the scale density, and the scale-length of the dark 
halo, respectively. The ratio of the virial radius (rvir) and 
Ts is a parameter in the NFW profile and described as c. A 
dark matter halo with the total mass of 6.5 x IO^Mq and 
c — 10 has Ts = 2.5 kpc and the virial velocity of 33 km 
s~^ in the NFW profile. The total mass in a NFW model is 
defined as the mass within Tvir in the present study. 



3.2 The Burkert profile 

The Burkert profile is described as: 

Pdm,0 



pdm(r) 



{r + ro){r^ + ro^) ' 



(2) 



where pdm,o and ro are the central dark matter density and 
the core (scale) radius, respectively. In the Burkert profile 
(B95), a relation between the total mass within the core 
(Mo) and ro is described as: 



^» = ^-2 X 10'( ig^)"''^- 



(3) 
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Figure 6. The same as Fig. 1 but for tlie model NFW3 in which 
the NEW profile with rs = 17.7 kpc and c = 5, J = 40°, and 
Mss/Ly = 2.3 are adopted. In this model, unusually large mass 
and low density of the dark matter halo for the SMC are assumed 
just for comparison with the results shown in Fig. 4 and 5. 



The virial radius (Rh) and the total mass of the dark matter 
halo in this profile are equal to 3.4 x ro and 5.8 x Mo, re- 
spectively. Therefore the total mass of the dark matter halo 
within Rh in the SMC can determine the core radius (ro). 
Fig. 3 shows a relation between ro and Mt (= Mdm + Mb) 
derived from the equation (3) for three different Ms/Ly. 
The parameter values (e.g., rs and ro) of models that are 
investigated in detail are summarized in the Table 1. 

The tidal radius of the SMC is about 5 kpc (Cardiner 
& Noguchi 1996) so that the shape of the dark matter halo 
can keep its (assumed) spherical shape within the central 3 
kpc after the LMC-SMC interaction. Therefore, the assumed 
spherically symmetric distribution of the dark matter halo 
in both the NFW and Burkert profiles can be regarded as 
reasonable even for the SMC about 0.2 Gyr after the last 
strong interaction with the LMC. 



3.3 Comparison with observations 



Fig. 4 shows that the model Bl having the Burkert profile 
with Mdm = 6.5 X IO^'Mq and ro = 3.2 kpc (i.e., Mt = 8.1 x 
lO^M© and R^ = 11.0 kpc) can well reproduce observations 
in terms of both (i) the slowly rising Vc profile in the inner 
part {R < 2 kpc) and (ii) the maximum T4 (~ 60 km s^^). 
The total mass of the dark matter halo within 3 kpc is about 
IO^Mq, which is still smaller than the total baryonic mass. 
These results imply that the SMC can have a large core in 
the mass distribution of the dark matter halo. We confirm 
that more massive models (e.g., B2) with larger ro can also 



reproduce well the observed K profile. The results of the 
models with lower Ms/Lv (=0.6 and 1) are shown in the 
Appendix A for comparison. 

Fig. 5 shows the simulated K as a function of R for the 
model NFWl having the NFW profile with c = 10, r^ = 2.5, 
and Mdm being exactly the same as that used in the Burker 
profile shown in Fig. 4. It is clear that although the NFW 
profile can reproduce the maximum K around i? ~ 3 kpc, it 
can not fit to the observed Vc profile as well as the Burkert 
one. In particular, in the central region {R < 2 kpc) the 
simulated profile is systematically well above the observed 
data points. The higher mass-density in the inner region 
of the dark matter halo described in the NFW profile is 
responsible for this discrepancy. These results imply that the 
NFW profile with reasonable model parameters are highly 
unlikely to reproduce the observations. 

It should be however stressed that au unusually low- 
density (i.e., smaller c), high-mass dark matter models with 
the NFW profile can reproduce the observations as well as 
the Burkert profiles. For example, the simulated Vc profile 
in the model NFW2 with Mdm = 6.5 x IO^Mq, c = 5, and 
rs = 5.1 kpc can be well consistent with the observations. 
Fig. 6 shows the simulated Vc profiles for the model NFW3 
with Mdm = 7.9 X 1O^°M0, c = 5, and rs = 17.7 kpc. In this 
extreme example with a low-density dark matter halo, the 
observations can be well reproduced, although it is unlikely 
that the mass of the SMC is much larger than that of the 
LMC (~ 2 X lO^^M©) that was previously considered to form 
a binary with the SMC (e.g., Murai & Fujimoto 1980). 

These results imply that the Burkert profile is better 
than the NFW one in terms of reproducing HI rotation curve 
in the SMC. Although these results do not rule out the NFW 
profile for the distribution of the dark matter halo of the 
SMC, the observed low-density dark matter halo (or large 
central concentration of the baryonic components) needs to 
be explained in the context of hierarchical clustering sce- 
narios. If, however, the SMC has such a large initial mass 
of Mdm ~ 7.9 X lO^^Mo, as shown in Fig. 6, then our views 
of the formation and evolution processes in the Magellanic 
system can be dramatically changed. 



4 DISCUSSION 

4.1 A massive past of the SMC 

If the SMC has low-mass and low-density of the dark mat- 
ter, as shown in our models, then the present results have 
the following four implications. Firstly, the SMC could have 
lost most of its original mass through tidal interactions with 
the Galaxy and the LMC owing to the extended initial dis- 
tribution of the dark matter halo. Essentially, the present- 
day SMC is a stripped core of an initially much larger 
galaxy. The SMC's tidal radius (rt) against the Galaxy can 
range from 5 kpc to 7.5 kpc for the total mass (Mt) of 
3 — 10 X IO^Mq (see the formula for the mass estimate in 
Gardiner & Noguchi 1995). Therefore i?h (H kpc) for the 
SMC with Mdm = 6.5 x lO^M© is significantly larger than 
rt, which means that the SMC could have lost all compo- 
nents outside rt owing to tidal stripping. The stripped mat- 
ters may well be now observed as the MS, the MB, and the 
stream of dark matter in the outer halo of the Galaxy. 
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Figure?. Orbital evolution of the MCs for Msmc = 8.1x10^"Mq 
(top), Msmc = 8.1 X 1O^M0 (middle), and Mamc = 3.0 X IO'^Mq 
(bottom). The time evolution of distances between MCs (thick 
solid), the LMC and the Galaxy (thin solid), and the SMC and 
the Galaxy (dotted) are separately shown. Here the total mass of 
the LMC is assumed to be 2.0 x IO^OM©. 



Secondly, the SMC could have much more strongly in- 
fluenced the LMC owing to its larger mass, in particular, 
for the last 3 Gyrs. The previous models assumed that 
Mt --- 3 X IO^Mq (e.g., Gardiner & Noguchi 1995; Yoshizawa 
& Noguchi 2003; Bekki & Chiba 2005, BC05), which thus 
would have underestimated gravitational influences of the 
SMC on the LMC. Thirdly, the SMC may well initially have 
a larger amount of HI gas owing to its large total mass. The 
total gas mass of the SMC in previous MS models is lim- 
ited to less than 10® M0 owing to the adopted models with 
smaller Mt ( e.g., Yoshizawa & Noguchi 2003). This limi- 
tation was responsible for the less consistent model for the 
total mass of the MS and the MB (Yoshizawa & Noguchi 
2003). Thus if future theoretical works adopt a more mas- 
sive SMC model, they will much more consistently explain 
the observed total mass (> lO^M©) of the MS and the MB. 

Fourthly, future orbital evolution models for the MCs 
would need to be revised significantly by considering the 
possible large total mass of the SMC. Since the orbital evo- 
lution of the MCs, which strongly depends on the mass of 
the SMC, is the most important ingredient for the forma- 
tion of the MS and the MB (e.g., Yoshizawa & Noguchi 2003; 
MuUer & Bekki 2007), the models with larger Mt can de- 
scribe the formation processes in a quite different manner: 
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Figure 8. The dependences of Mo/Ly on My (upper) and on 
a (lower) within effective radii for the dwarfs in the Local group 
(open squares) and the SMC (filled circles). Here Me and a repre- 
sent the total masses (including dark matter halos) within effec- 
tive radii and velocity dispersions for dwarfs, respectively. Here 
the observational results by Mateo (1997) are used and the solid 
line in the lower panel is from Zaritsky et al. (2006). 



future theoretical works on the MS and MB formation need 
to investigate models with a larger initial mass of the SMC. 

In order to illustrate how the orbital evolution of the 
MCs can be changed owing to a larger mass of the SMC, we 
investigate the evolution based on the same orbital evolu- 
tion model used in previous works (e.g., BC05): the adopted 
initial parameters on initial positions and velocities of the 
MCs and the gravitational potential of the Galaxy are ex- 
actly the same as those in BC05. Fig. 7 shows the orbital 
evolution of the MCs for the models in which (i) the meth- 
ods to model orbital evolution of the MCs are the same as 
those used in BC05 as "no dynamical friction" , and (ii) Mt 
{— Msmc) is equal to or much larger than that used in BC05 
(3 X 10® M0). It is clear from this figure that even if dy- 
namical friction between the MCs is not included, the MCs 
can not keep their binary status for more than a few Gyr in 
the model with large Mt- We confirm that if the SMC has 
Mt > 1.5 X lO^^M©, the MC can not keep its binary status 
for more than 4 Gyrs. 

If dynamical friction between the MCs can be effective, 
the MCs can not keep their binary status for more than 4 
Gyrs even for Mt = 3 x 10^°Mq (BC05): the more massive 
SMC suggested in the present study is highly unlikely to be 
coupled with the LMC for more than 4 Gyrs. Furthermore, 
the evolution processes (e.g., star formation history) of the 
SMC can be significantly different in the models with larger 
Mt owing to the stronger self-gravity of the SMC against 
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tidal perturbations from the Galaxy and the LMC. Thus 
HI observations on the mass distribution of the SMC (e.g., 
S04) may well change significantly views about formation 
and evolution of the MCs. 

Sofue (1998) analyzed the rotation curve profile of M82 
and found that the profile declines in a Keplerian fash- 
ion outside 200 pc. He thus suggested that (i) a significant 
amount of mass is missing in the outer part of M82 and (ii) 
tidal stripping due to past tidal interaction between M81 
and M82 is responsible for the missing mass. These results 
by Sofue (1998) and the present ones suggest that the SMC 
would be the second case of low-mass galaxies that have lost 
a significant fraction of their original masses. A significant 
difference between these two galaxies is that the rotation 
curve of the SMC is not observed to show a decline in a 
Keplerian fashion even for the outer part. 



nents within the central 3kpc, then the origin of this unique 
nature of the SMC needs to be understood. Efficient trans- 
fer of gas to the central regions of galaxies and the resultant 
starbursts there in galaxy merging (e.g., Mihos & Hernquist 
1995) may significantly increase the baryonic fractions of 
merger remnants. Recently, Bekki & Chiba (2008) have sug- 
gested that the SMC was formed from merging of two gas- 
rich dirr's with extended HI gas very long time ago. There- 
fore, the observed large fraction of baryonic component in 
the SMC can be due mainly to the past major merger event 
that formed the SMC. Thus the observed HI kinematics of 
the SMC can provide unique insights into the formation his- 
tory of the SMC. 



4.2 The closest dirr or dSph with a flat DM core 

The present study has first suggested that the dark mat- 
ter halo of the SMC has a flat core well described by the 
Burkert profile rather than the NFW one: the SMC can be 
the nearest example of dirr's (or dSph's) in which the NFW 
profiles are hard to fit to HI rotation curves. It should be 
however stressed that Valenzuela et al. (2007) have recently 
suggested a possibility of a large underestimate of the HI 
rotational velocity in the central 1 kpc region of two dwarf 
galaxies due to noncircular motions. Since Valenzuela et al. 
(2007) adopted rotating disk models for their investigation, 
their results can not be applied directly to the SMC which 
has no evidence of stellar rotation (Harris & Zaritsky 2006) . 

However, if the last close LMC-SMC interaction (Bekki 
& Chiba 2007) could cause the radial flow of HI gas and 
thus the noncircular motions, then the observed K in S04 
could be underestimated. If this is the case, the NFW pro- 
file with a cuspy core would be still possible. The present 
models are not based on numerical simulations of HI gas in 
the SMC thus can not allow us to make it clear whether the 
underestimation of circular velocities suggested by Valen- 
zuela et al. (2007) is equally possible in the gas disk of the 
SMC interacting with the Galaxy and the LMC. We plan 
to discuss this problem in our future paper in which the 
HI rotation curve of the SMC will be constructed based on 
hydrodynamical simulations of the SMC evolution. 

The derived small mass of the dark matter halo within 
the optical radius of the SMC suggests that the SMC is a 
baryon-dominated galaxy. This strong central concentration 
of baryonic components in the SMC also suggests that the 
formation history of the SMC can be significantly different 
from that of other dirr's dominated by dark matter halos. 
Fig. 8 shows how unique the SMC is by comparing the lo- 
cation of the SMC and dwarfs in the Local group on the 
Mv — Mc/I/v plane and a — Mc/Lv one, where Me and 
a represent the total masses (including dark matter halos) 
within effective radii and velocity dispersions, respectively. 
Although the SMC has much lower Mv/Mc in comparison 
with other dwarfs in the Local group, it is possible that 
My /Ma is not so low for a luminous dwarf population (i.e., 
more luminous dwarfs can have lower My/Mc; see Zaritsky 
et al. 2006). 

If the SMC really has a very low-density dark matter 
halo and thus a strong concentration of its baryonic compo- 



4.3 Stronger constraints on the mass of the SMC 
from future observations 

One of potential problems in estimating the mass of the 
SMC based on the observed Vc profile is that we can un- 
derestimate significantly the SMC's inclination angle (e.g., 
by 20 degrees in i). The HI distribution and kinematics are 
highly disturbed by numerous expanding shells (S04), mak- 
ing a more accurate estimate of the inclination angle very 
difficult. A significant underestimate of the inclination an- 
gle, although highly unlikely, would result in a larger total 
mass of the SMC within the central 3 kpc (e.g., 5.3 x lO^M© 
for the true 14 of 80 km s~^). Therefore, future observations 
of the distribution and radial velocities of a large number of 
older stars (e.g., RGB and AGB stars) are very important 
to derive the mass of the SMC independently from the HI 
observations (S04). 

Although the present study has shown that the SMC 
can have a small amount of dark matter within the central 
3 kpc, the detailed profile of the dark matter distribution in 
the SMC has not been determined. The observed Vc profile 
within the central 3 kpc alone does not allow us to distin- 
guish between the Burkert model with Mt --^ 6.5 x lO^M© 
and a very massive Burkert (or even NFW) profile with 
Aft > 10^°Mq. The mass profile of the SMC at i? > 3 kpc, 
where less massive models can show almost flat Vc profiles, 
need to be observationally investigated to provide stronger 
constraints on the distribution of the dark matter halo in 
the SMC. 

Noel & Gallart (2007) have recently found intermediate- 
age and old stars belonging to the SMC but being located 6.5 
kpc from the SMC center. This suggests that the SMC stellar 
population may be more massive than previously thought. 
Future spectroscopic observations of radial velocities of the 
stars in the outer halo will enable us to derive the mass pro- 
file for the outer regions of the SMC and thus to better con- 
straint the radial mass profile of the dark matter halo. The 
ongoing VMC project (the VLT VISTA near-infrared YJKs 
survey for the MCs), combined with other large projects, 
will soon provide structural and kinematical data sets for 
the outer regions of the SMC (e.g., Cioni et al. 2007). We 
plan to discuss the entire mass profile of the dark matter 
halo of the SMC based both on these ongoing observations 
and previous HI observations (S04) in our future papers. 



The lost dark matter 9 



5 CONCLUSIONS 

We have investigated the mass (A/dm) and the mean den- 
sity (pdm) of the dark matter halo within the optical radius 
(~ 3 kpc) of the SMC based on the high-resolution HI ob- 
servations of the SMC. Although there are some limitations 
of the present models in estimating the total mass of the 
SMC solely from the observed HI rotation curve by S04, we 
have tried to derive the mass and best possible dark matter 
models. The results and their implications are described as 
follows. 

(1) The total mass, mean density, and distribution 
of the dark matter halo depend strongly on the adopted 
AIs/Ly. For example, the estimated Mdm and pdm are 
7.5 X IO^'Mq and 6.7 x lO'^M© pc~^, respectively, for 
Ma/Lv = 3.8. This implies that if the SMC is dominated by 
older stellar populations, the SMC has a low-density dark 
matter halo described by B95 rather than by the NFW. 

(2) The models with lower Mg/Ly (< 1) can show 
higher Afdm and pdm yet have difficulties in reproducing the 
observed rotation curve, in particular, the inner part {R < 1 
kpc) of the curve, both for the Burkert and the NFW pro- 
files. This implies that the SMC is unlikely to be dominated 
by younger stellar populations. 

(3) The Burkert profile can fit better the observed HI 
rotation curve of the SMC than the NFW profile for a rea- 
sonable Afdm and Afs/Ly for reasonable Ms/Ly (~ 2). How- 
ever the NFW profiles with low-densities, small c (~ 5), and 
large total masses (Afdm > 10^°Mq) can be also fit well 
the observations, though such low-density NFW profiles are 
highly unhkely for low-mass dwarfs hke the SMC. The SMC 
can be the closest example among dirr's in which rotation 
curves are difficult to fit to the mass profiles of dark matter 
halos predicted by hierarchical clustering models. 

(4) The dark matter halo of the SMC is likely to have 
the initial total mass of Mdm ~ 6.5 x lO^M© and the core 
radius of ro ~ 3.2 kpc, which means that the SMC could be 
much more massive than what was previously thought. 

(5) Orbital evolution models of the SMC and the LMC, 
on which formation models of the MS and the MB are based, 
would need to be revised by considering the derived much 
larger total mass of the SMC. The observed large concen- 
tration of the baryonic components in the SMC is due to its 
unique formation history and is very different from those of 
other dirr's. 
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APPENDIX A: THE MODELS WITH LOWER 
MASS-TO-LIGHT-RATIOS 

The models with lower Ms/Ly (i.e., smaller stellar masses) 
have larger masses of the dark matter halo and thus larger 
core (or scale) radii for the SMC. It would be therefore in- 
teresting to discuss whether such more massive models can 
explain the mner part of the observed rotation curve. Fig. 
Al shows that for M^/Ly = 0.6, even the very massive 
dark matter model (B3) with ro = 9.5 kpc (corresponding 
to A/dm = 7.9 X 10^°Mq) with the Burkert profile can not 
explain the observed profile of the rotation curve owing to 
the required large mass in the inner region of the SMC. Fig. 
A2 also shows that for M^/Ly = 1.0, the massive model 
(B4) can not reproduce the inner part of the observed ro- 
tation curve owing to the low density of the halo. It should 
be stressed here that the SMC is highly unlikely to have a 
mass much larger than that of the LMC (~ 2 x 10^°Mq. 

These results imply that in order for the observed rota- 
tion curve to be reproduced by the present model, a higher 
Ms/Ly needs to be adopted: the SMC is highly unlikely to 
be dominated by younger stellar populations. This is con- 
sistent with the results by Harris & Zaritsky (2004) who 



10 K. Bekki and S. Stanimirovic 




R (kpc) 

Figure Al. The same as Fig. 1 but for the model B3 in which 
the Burkert profile with ro = 9.5 kpc, i = 40°, and Ms/Ly = 0.6 
are adopted (M^m = 7.9 X IO^OMq). 




R (kpc) 

Figure A2. The same as Fig. 1 but for the model B4 in which 
the Burkert profile with ro = 9.5 kpc, i = 40°, and Ms/Ly = 1 
are adopted (Mdm = 7.9 X IO^OMo). 



showed that about 50% of the stellar populations in the SMC 
were formed prior to 8.4 Gyr ago. Given that the large frac- 
tion of stars in the SMC are suggested to be formed between 
2.5 and 3 Gyr ago in starbursts (Harris & Zaritsky 2004), 
the SMC must be dominated by intermediate-age and old 
stellar populations. Thus we can conclude that our success- 
ful reproduction of observations in the models with higher 
Ms/Ly ~ 2 — 4 is quite reasonable. 



